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Abstract 


We present here a newly developed satellite cloud retrieval algorithm SACURA. The algorithm is based on the modified expo- 
nential approximation for the cloud reflection function. It is valid for clouds having the optical thickness larger than 5. Retrieval 
results for selected cloudy fields using data from SCIAMACHY onboard ENVISAT are presented. The cloud algorithm itself does 
not use the look-up-table approach and can be applied to other instruments placed on satellite and airborne platforms. 
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1. Introduction 


Clouds play an important role in the Earth climate 
system (Kondratyev and Binenko, 1984; Liou, 1992). 
The amount of radiation reflected by the Earth—-atmo- 
sphere system into outer space depends not only on 
the cloud cover and the total amount of condensed 
water in the Earth atmosphere but also on the size of 
droplets der, the thermodynamic state of water in clouds, 
and vertical distributions of cloud parameters. The 
information about microphysical properties and spatial 
distributions of terrestrial clouds on a global scale can 
be obtained only with satellite remote sensing systems. 
Different spectrometers and radiometers (Bovensmann 
et al., 1999, 2004; Deschamps et al., 1994; King et al., 
1992; Nakajima et al., 1998), deployed on space-based 
platforms, measure the angular and spectral distribution 
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of intensity and polarization of backscattered solar light. 
Generally, the measured values depend on both geomet- 
rical and microphysical characteristics of clouds. Thus, 
the inherent properties of clouds can be retrieved (at 
least in principle) solving the inverse problem. The accu- 
racy of the retrieved values depends on the accuracy of 
measurements and forward radiative transfer model. 
In particular, it is often assumed that clouds can be rep- 
resented as homogeneous and infinitely extended in the 
horizontal direction plane-parallel slabs (Goloub et al., 
2000; Han et al., 1994; King, 1987; Nakajima and King, 
1990; Nakajima et al., 1991; Rossow, 1989; Rossow and 
Schiffer, 1999; Kokhanovsky et al., 2003). The range of 
applicability of such an assumption for real clouds is 
very limited as is shown by observations of light from 
the sky on a cloudy day. For example, it was found that 
the retrieved cloud optical thickness t depends on the 
viewing geometry (Loeb and Davies, 1996; Loeb and 
Coakley, 1998). This, of course, would not be the case 
for an idealized plane-parallel cloud layer. However, 
both the state-of-art radiative transfer theory and com- 
puter technology are not capable to incorporate three 
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dimensional effects into operational satellite retrieval 
schemes. As a result, cloud parameters retrieved should 
be considered as a rather coarse approximation to real- 
ity. However, even such limited tools produce valuable 
information on terrestrial clouds properties. For exam- 
ple, it was confirmed by satellite measurements that 
droplets in clouds over oceans are usually larger than 
those over land (Han et al., 1994). This feature, for in- 
stance, is of importance for the simulation of the Earth’s 
climate (Slingo, 1989). Microphysical characteristics of 
clouds are of importance for dynamical meteorology. 
For example, the probability of precipitation increases 
for larger effective radii. 

A new Semi-Analytical CloUd Retrieval Algorithm 
(SACURA) for the cloud liquid water path (LWP) 
and water droplets size determination presented here is 
based on the asymptotical solution of the radiative 
transfer equation for a special case of disperse media, 
having large optical thickness. This solution was ob- 
tained by Germogenova (1961, 1963) for plane-parallel 
turbid slabs. Such an approach has already been used 
in a number of studies (Rozenberg et al., 1978; King, 
1987). The difference is that the asymptotical solutions 
are further simplified such that the inverse problem is re- 
duced to the solution of a single transcendent equation. 
This allows us to speed up the retrieval process signifi- 
cantly without the substantial loss of the accuracy of 
the retrieved cloud characteristics. 

The algorithm is restricted to the case of optically 
thick clouds (the optical thickness t > 5). It is planned 
to be supplemented in future by the look-up-table 
(LUT) approach at t<5 (Nakajima and King, 
1990). However, it should be stressed that the optical 
thickness is highly correlated with the geometrical 
thickness of clouds. For instance, it was shown (Fei- 
gelson, 1981) that clouds having t<5 have the geo- 
metrical thickness less than 200m on average. It is 
difficult to expect that such clouds are homogeneous 
in the horizontal direction. Thus, the horizontal pho- 
ton transport must be accounted for (Platnick, 2001). 
However, this is not considered in standard retrieval 
procedures (Arking and Childs, 1985; Nakajima and 
King, 1990; Platnick et al., 2001, 2003; Nakajima 
et al., 2005). So results obtained from more accurate 
forward modelling as compared to ours are anyway 
biased due to cloud inhomogeneity effects. This gives 
another justification besides speed concerns to the ap- 
proach we use. 

Besides dep and t, we also retrieve the liquid water 
path w (gm~’), the cloud albedo r, and the column con- 
centration of droplets N (m~°), using well-known rela- 
tionships between these quantities (see below). 

The cloud top height is derived from measurements in 
the oxygen A-absorption band as recommended by 
Yamomoto and Wark (1961). Such an approach has 
been used extensively by many authors (Kuze and 


Chance, 1994; Fischer et al., 2000; Koelemeijer et al., 
2001, 2002; Kuji and Nakajima, 2002). 

Finally, the cloud thermodynamic state is obtained 
using the fact that spectral signatures of liquid water dif- 
fer from those for ice in the selected spectral range 1550- 
1670 nm (Pilewskie and Twomey, 1987a,b; Knap et al., 
2002; Acarreta et al., 2004). 


2. The determination of the cloud optical thickness, the 
cloud albedo, the liquid water path, the column droplet 
concentration, and the effective size of particles using 
backscattered solar light measurements 


The SACURA is based on the two-wavelength semi- 
analytical cloud retrieval algorithm for the liquid water 
path w and the cloud particle size ag determination pro- 
posed by Kokhanovsky et al. (2003). It is known that 
the reflection function in visible is primarily dependent 
on the cloud optical thickness (and, therefore, on the 
LWP) and the infrared channels are sensitive to the size 
of particles in clouds (Nakajima and King, 1990). 

Parameters dep and w are used to find the cloud opti- 
cal thickness t, the cloud albedo r, and the column drop- 
let concentration N (cm~’) from following approximate 
analytical equations valid in the visible (Kokhanovsky, 
2004a,b): 


1.5w + 1.1 1 
T= r= 
paet (kaaf J’ 1.07 + 0.751(1 — g)’ 
E 
~ 2nta2,’ 


(1) 


where p = 1 gcm_° is the density of water, k = 27/4, A is 
the wavelength, ¢ = 0.7 and 


g = 0.08 — 0.5(kaer) °’, (2) 
is the asymmetry parameter. 
The values of agp and w are found from the system of 


two algebraic equations (Kokhanovsky et al., 2003; 
Kokhanovsky, 2004c): 


ti (der, w)[1 — Aj] 


Ri =R? 
i © 1—A[l — t (as, w)] 


Ko(¥)Ko(¥o), (3) 


Ro = RY, exp(—V2(aer)(1 — 0.05y2(aer))u) 
ty (Ger, w) 
1 = Aory( Ger, w) 
X th(der, W)Ko(V)Ko(Vo)- (4) 


— |exp(—x2 (der, w) — ¥2(Ger)) 


The subscripts “1” and “2” refer to wavelengths 4; and 
22 in visible and near - infrared channels, respectively. 
R, and R, are measured reflection functions in the visi- 
ble and infrared, respectively. Specifically, the SACURA 
uses measurements at the wavelengths 443 and 1550 nm. 
These channels are almost free of gaseous absorption. 


A.A. Kokhanovsky et al. | Advances in Space Research 36 (2005) 789-799 791 


Table 1 
Auxiliary functions (Kokhanovsky and Rozanov, 2003; Kokhanovsky et al., 2003; Kokhanovsky, 2004a; Van de Hulst, 1980) 
The auxiliary function Formula 
u Ko(V)Ko(vo) 
R? (8, vo, p) 


Ko) 
R? 
ty 


r2 
X2 
y2 


y2 


3 (1 + 2cosð) 
3.944 — 2.5(cos Vo + cos 9) + 10.664 cos 9o cos 0 + F (0) 
4(cos 0 + cos Jo) 
sinh y, 4.86 — 13.08 cos vo cos V + 12.76 cos W cos 7 
sinh[x2 + 1.07y,] T 
exp(=y2(1 — 0.05y2)) — t2 exp(—x2 — ya) 
Yat 
4y2/3(1—g2) 
3(1 — g2)(1 — %2) 


exp[xa] 


Here, R? is the reflection function of a semiinfinite nonabsorbing cloud layer, Ko(1) is the escape function, which describes the angular distribution of 
light transferred by a thick cloud layer, t; and r; are hemispherical cloud transmittance and reflectance, respectively, for diffuse light illumination 
conditions, and y; is the diffusion exponent. The index j= 1, 2 corresponds to the spectral channel used. The function F(@) is calculated as 
F(0) = Sop) Xm[Pm(cos 0) — (— 1)” Pm(cos 7)Pm(cos vo)], where Xm are coefficients of the expansion of the phase function with respect to Legendre 
polynomials P,,,(cos 0) (Kokhanovsky, 2004b). Here 0 = cos~!(—cos Vecos Vo + sin J sin Vocos o) is the scattering angle, Vo is the solar zenith angle, 9 is 
the observation zenith angle, and ¢ is the relative azimuth. The value of ft, can be calculated from r given by Eq. (1) at the wavelength 22. Namely, it 


follows: t; = 1 — r(4,). 


Values of A, and A> give us the surface albedos in the 
visible and near-infrared, respectively, Vo is the solar ze- 
nith angle, V is the zenith observation angle. Other func- 
tions in Eqs. (3) and (4) are given in Table 1. 

We underline that the land surface reflectance gener- 
ally decreases for shorter wavelengths. This allows to ne- 
glect its influence on retrievals for optically thick clouds 
considered in this paper at the wavelength 443 nm. Our 
simulations show that for bright surfaces like snow and 
ice the surface reflectance must be accounted for in the 
retrieval procedures even for relatively optically thick 
clouds having values of t below 70 for water clouds 
and 40 for ice clouds. 

Eqs. (3) and (4) have two unknowns (aef and w). 
Standard methods and programs are available for the 
solution of such problems. However, taking into ac- 
count the specific form of Eq. (3), the exclusion method 
is used here. So we can formulate a single transcendent 
equation with one unknown (aer). It follows from Eq. 
(3) that 


n = (ERK A I 
R? — R; 1-A,Jf ` 


(5) 


On the other hand, it follows from the formula for tı gi- 
ven in Table 1 in the visible: 


4(t;! — 1.07) 
3(1 — g(aer)) ’ 
where we neglect A (see Table 1) for the sake of simplic- 
ity and the parameter g, is given by Eq. (2) at 2 = 2. 
The optical thickness q; is related to the liquid water 


path w and the effective radius of droplets by the follow- 
ing equation (see Eq. (1)): 


T = WGext(A1, aef), (7) 


Ti = 


(6) 


where 


L5 Ll 
Gea = —— [1+ . 8 
l paa ( (kacr)? ) % 


Therefore, the analytical expression for w can be ob- 
tained from Eqs. (6)-(8): 
—1 
pa 4(¢;! — 1.07) (9) 
3(1 — & (er))Gext(A1, Ger) 

where the dependence of all functions on the value of the 
effective radius is explicitly noted. Thus, one unknown 
parameter (w) is expressed in terms of another (aer). 
The substitution of Eq. (9) into Eq. (4) with account 
for formulae in Table 1 yields a single transcendent 
equation for the determination of the effective radius 
of droplets in a cloud. It should be pointed out that w 
enters Eq. (4) via X2(dep,w). The parameter 


3(1-g2) 
nition (see Table 1). Here, wo2 and gə are the single scat- 
tering albedo and the asymmetry parameter at À = A), 
respectively. The function x2(dep, w) can be written as : 


x(der, w) = V/3(1 — @o2(aer))(1 = 85(der))t2 (Ger, w), 
(10) 


where only t2 depends on the liquid water path w. The 
index “2” corresponds to the wavelength A, in the 


y, = 4, / 722 in Eq. (4) does not depend on w by defi- 


near-infrared. We use following results valid at 
Az = 1550 nm : 
4 n 
1= g = X a(2rag 4)’, 
i=0 
= _ azn hy a1) 
Sa 22pGext(A2, def) - cia A2) 
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Table 2 

Parameters c; and d; in Eq. (11) 

i 0 1 2 3 4 

Gi 0.0608 2.465 —32.98 248.94 —636 


d; 1.671 0.0025 —0.0002365 2.86e—6 


where parameters c; and d; are given in Table 2. These 
expressions were obtained fitting Mie results for spheri- 
cal polydispersions with account for general geometrical 
optics results for large spherical droplets (Kokhanovsky, 
2004a). 

Taking into account the definition of the optical 
thickness: t = dextL, where L is the cloud geometrical 
thickness and Gext is the extinction coefficient, t2 in Eq. 
(10) is given by: t2 = ®(ae¢)t,. The analytical expression 
for the function 
Sext (Az) 
Plar) Gext(A1) f 
is easily derived from Eq. (8). 

Summing up, the substitution of Eq. (10) and also 
other parameters from Table | into Eq. (4), results in 
a single equation for the effective radius retrieval from 
measurements of R; and R2. This equation can be easily 
solved numerically (Press et al., 1992). The obtained va- 
lue of der is used to calculate the optical thickness via Eq. 
(1), where the liquid water path is determined directly 
from Eq. (9) with account for Eq. (5). The values r 
and N are found using Eq. (1). Interestingly, it is possi- 
ble to obtain the spherical albedo r from bi-directional 
reflectance measurements. 


(12) 


3. The cloud top height determination 


The determination of the cloud top height using 
SACURA is based on measurements of the top-of- 
atmosphere (TOA) reflection function R in the oxygen 
A-band (Kokhanovsky et al., 2004, 2005a,b; Rozanov 
and Kokhanovsky, 2004; Rozanov et al., 2004). The 
cloud reflection function is extremely sensitive to the 
cloud top height in the center of the oxygen absorption 
A-band (Yamomoto and Wark, 1961). To find the value 
of h, we first assume that the TOA reflectance R can be 
presented in the form of a Taylor expansion around the 
assumed value of the cloud top height equal to ho: 


RA) = Rho) + Y a-h, (13) 


where a; = R(ho)/i! Here, R(Ao) is the i-derivative of R 
at the point ho. The next step is the linearization, which 
is the standard technique in the inversion procedures 
(Rozanov et al., 1998; Rozanov and Kokhanovsky, 
2004; Rozanov et al., 2004; Rodgers, 2000). We found 


that the function R(/) is close to a linear one in a broad 
interval of the argument change (Kokhanovsky and 
Rozanov, 2004). Therefore, we neglect nonlinear terms 
in Eq. (13). Then it follows: 


R = R(ho) + RB'(ho) (A — ho), (14) 


where R’ = R, We assume that R is measured at several 
wavelengths in the oxygen A-band. Then instead of the 
scalar quantity R we can introduce the vector Ries With 
components (R(21), R(A2),...,RUn)). The same applies 
to other scalars in Eq. (14). 

Therefore, Eq. (14) can be written in the following 
vector form: 


y=xd (15) 


where y = Recs R(ho), a= R (ho), and x = h — ho. Note 
that both measurement and model errors are contained 
in Eq. (15). The solution £ of the inverse problem is ob- 
tained by the minimizing the following cost function 
(Rodgers, 2000): 


® = ||y — xa’, (16) 


where || || means the norm in the Euclid space of the cor- 
respondent dimension. 

The value of x, where the function ® has a minimum 
can be presented as (see Eq. (15)): 


(17) 


where (,) denotes a scalar product in the Euclid space. 

Therefore, from known values of the measured spec- 
tral reflection function Rmes (and also values of the calcu- 
lated reflection function R and its derivative R’ at h = ho) 
at several wavelengths, the value of the cloud top height 
can be found from Eq. (17) and the equality: h = £ + ho. 
The value of ho can be taken equal to 1.0 km. This is a 
typical value for low level clouds (Feigelson, 1981). The 
main assumption in our derivation is that the depen- 
dence of R on h can be presented by a linear function 
on the interval x (Kokhanovsky and Rozanov, 2004). 

As a matter of fact, the SACURA code finds both the 
cloud top height h and the cloud geometrical thickness / 
simultaneously. This requires the minimization of the 
following cost function (see Eq. (16)): 


® = ||y — AX)’. (18) 


The elements of the matrix 4 are correspondent weight- 
ing functions (Rozanov et al., 1998, 2002; Rodgers, 
2000). The solution of the inverse problem is given by 
the vector-parameter ¥. This vector has five compo- 
nents, which give corrections to the initially assumed 
cloud top height and cloud geometrical thickness, the 
correction to the initially assumed half-width of the 
spectrometer spectral response function, the shift 
parameter, and the squeeze parameter. Further details 
are given by Rozanov and Kokhanovsky (2004). 
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We have developed two versions of the retrieval algo- 
rithm. One is based on the exact radiative transfer calcu- 
lations of the reflection function R and its derivatives R’ 
and another one is based on the approximate represen- 
tation of R by the following equation (Kokhanovsky 
and Rozanov, 2004): 


R = Ro + TIRT, (19) 


where Ro gives the reflection function of the atmospheric 
layer above the cloud and R, is the cloud reflection. 
Functions T; (i=1, 2) give transmission coefficients 
from the sun to a cloud and from the cloud to a satellite, 
respectively. Approximate analytical equations for all 
functions in Eq. (19) are given by Kokhanovsky and 
Rozanov (2004). We have used the correlated k-distribu- 
tion method to account for the high-frequency oscilla- 
tions of the oxygen molecule absorption cross-section 
Ca (Lacis and Oinas, 1991; Buchwitz, 2000; Buchwitz 
et al., 2000). The temperature and pressure dependence 
of o, for a given location of measurements was ac- 
counted for using the standard atmosphere model built 
in SCIATRAN (Rozanov et al., 2002). We have used 
the most recent version of the HITRAN molecular spec- 
troscopic database to get data on cross-sections o, 
(Rothman et al., 2003). 

For the operational retrievals, we have used the re- 
trieval technique based on the approximate representa- 
tion of R (see Eq. (17)). It allows to speed up the 
retrieval process considerably. We underline that the 
SACURA is capable to retrieve the cloud geometrical 
thickness as well. However, the meaning of this param- 
eter for multi-layered cloud systems is not clear. 


4. The phase index 


The thermodynamic state of clouds is an important 
parameter for cloud physics, climate problems and avi- 
ation safety. It is known that spectral properties of ice 
differ significantly in the spectral range 1500-1700 nm 
from those of water (Knap et al., 2002; Acarreta et al., 
2004). Therefore, we introduce the phase index (PI) «: 


_ R(1550 nm) 


* = R(1670 nm)’ (20) 


to characterize the thermodynamic state of a cloud. Low 
values of this parameter correspond to ice clouds. 

We find using radiative transfer calculations that the 
value of « is in the range 0.7-1.0 for water clouds and it 
is in the range 0.5-0.7 for ice clouds. Calculations were 
performed using the modified asymptotic equations gi- 
ven by Kokhanovsky and Rozanov (2003) in the 
assumptions that both crystals and droplets have the 
spherical shape. It was assumed that the cloud optical 
thickness is varied in the range 5-30 for both water 
and ice clouds. The effective radius of droplets in water 


clouds and also ice crystals was changed in the range 
5-30 microns. Studies performed by Kokhanovsky 
et al. (2005b) show that the increase of der leads to the 
increase of the separation of both cloud types. Note that 
ice crystals have highly irregular shape and also the size 
is often in the range 100-500 microns as compared to 
just 10 microns for warm water clouds. This will en- 
hance differences even further. Therefore, the PI indeed 
can be used as an indicator of the cloud phase. More dis- 
cussions on the phase index are given by Kokhanovsky 
et al. (2005b). 

Obviously, we have for clear sky over black sur- 
faces: R(1550 nm) > R(1670 nm), and, therefore, « > 1. 
Some highly reflecting soils could have values of 
æ = 0.5-0.7 similar to crystalline clouds. However, these 
pixels are screened out using information on the de- 
rived cloud top height, the cloud optical thickness or 
both. 


5. Results of retrievals 


The SACURA has been applied to all nadir data of the 
SCAnning Imaging Absorption spectrometer for Atmo- 
spheric ChartographY (SCIAMACHY) on board ENVI- 
SAT. SCIAMACHY is a moderate resolution imaging 
spectrometer. It composed of a mirror system, a telescope, 
a spectrometer, and thermal and electronic sub 
systems. Spectra are recorder simultaneously from 214 
to 1750nm and also in two smaller spectral regions 
1940-2040 and 2265-2380 nm. Spectral resolution 
depends on the channel: 0.21 nm (214-314 nm), 0.22 nm 
(309-4004 nm), 0.47nm (392-605nm), 0.42nm 
(598-790 nm), 0.55 nm (776-1056 nm), 1.56 nm (991- 
1750 nm), 0.21 nm (1940-2040 nm), 0.24nm (2260- 
2384 nm). Measurements of SCIAMACHY are 
performed in nadir, limb, and solar/lunar occultation 
modes (14 orbits per day). The ground pixel size of 
SCIAMACHY varies depending on the spectral channel 
used. We used pixels with the size 30*60 km? in the retri- 
evals presented in this paper. Further details related to 
SCIAMACHY are given by Bovensmann et al. (1999). 
The results of retrievals and correspondent files are freely 
available at www.iup.physik.uni-bremen.de/scia-arc. 

The results of cloud top height retrievals using the 
oxygen A-band depth technique are given in Fig. 1 for 
September 17th, 2003. All 14 SCIAMACHY orbits 
available for this date are shown. The enhanced picture 
of every state is given at our website. The instrument 
operates in the limb mode in the case of gaps shown 
in the figure. Then nadir data are not available. White 
colour corresponds to cloud top heights above 15 km 
and the situation with no colours in the SCIAMACHY 
ground scenes corresponds to a clear sky conditions. It 
follows from Fig. 1 that the number of clear sky situa- 
tions is very limited in comparison with cloudy scenes. 
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Fig. 1. The global cloud top height map for September 17th, 2003 as derived by SACURA using SCIAMACHY data. 


This really should be the case because cloudiness covers 
approximately 60% of the terrestrial surface at any given 
time. So cloud retrieval algorithms are crucial for 
SCIAMACHY trace gases columns retrievals. 

It follows from Fig. 1 that most of high-altitude 
clouds are located in the region 0—30°N degrees. This 
coincides with general knowledge on global cloud verti- 
cal distribution available at present. Also hurricane 
activity is most pronounced namely in this region. In 
particular SCIAMACHY captured the Hurricane Isabel 
in the scene shown in Fig. 1. The Isabel occupies a part 
of the scene located North-East of Florida (30°N, 70°E) 
for the SCIAMACHY orbit 08095. A small portion of 
the cloud system associated with Isabel is also seen in 
the next state (in direction to the South). The distribu- 
tion of the cloud top heights with the latitude for the or- 
bit 08095 of ENVISAT, which contains data for Isabel, 


is given in Fig. 2. It follows that Isabel reaches height 
18 km. This is quite unusual for hurricanes at such lati- 
tudes to penetrate to such heights. 

To confirm our finding we analysed SCIAMACHY 
limb state just before the hurricane. Then SCIAM- 
ACHY has performed a vertical scan of the hurricane 
approaching it from the north. In particular, we plot 
the colour ratio c= N(1050 nm)/N(750 nm), where 
N(A) is the backscattered signal at the correspondent 
wavelength A normalized at its background value at 
36 km, in Fig. 3. It follows that this ratio is constant un- 
til the height 18 km. Then it drops substantially. It can 
be explained by the fact that both ice crystals and water 
droplets scatter light almost independently at the wave- 
length in the spectral interval 400-1200 nm. This is not 
the case for a clear atmosphere above the hurricane. 
Then the value of c decreases (e.g., c should be lower 


A.A. Kokhanovsky et al. | Advances in Space Research 36 (2005) 789-799 795 


20 


18 


16 


2k = = 
fo) N > 


œ 


cloud top height, km 


-80 -60 -40 -20 0 20 40 60 80 
latitude, degrees 


Fig. 2. The dependence of the cloud top height on the latitude for the 
orbit 08095 (September 17th, 2003) as derived by SACURA. 


for both stratospheric aerosol and molecular scattering 
than for ice or water clouds). This is seen in experimen- 
tal data presented in Fig. 3 as well. Our considerations 
confirm that the hurricane under study penetrates to a 
high altitude and signifies an important role of hurri- 
canes in supplying water and also other lower atmo- 
sphere constitutes into stratosphere. Hurricanes could 
also alter the stratospheric dynamics (at least in a hurri- 
cane season). 

Note that top of a hurricane should consist of ice 
crystals. This is confirmed by the phase index values ob- 
tained in this case. They are generally below 0.7. The 
phase index statistical distribution for a whole orbit 
08095 is shown in Fig. 4(a). It follows that the phase in- 
dex is in the range 0.35-1.0. It has a bi-modal distribu- 
tion. The mode positioned at values of PI «~ 0.85 is 
due to water clouds, and yet another mode located at 
the phase index close to 0.7 is due to crystalline clouds. 
The visual inspection of Fig. 4(a) confirms that the prob- 
ability of ice cloud occurrence is at least 2 times smaller 
than that of water clouds for the orbit analysed, which is 
a reasonable estimate. The statistical distributions of 
cloud top heights for the same orbit is shown in Fig. 
4(b). It follows that this distribution is bi-modal with 
most frequent values of the cloud top height located at 
5 and 7 km. The statistical distribution of the derived 
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Fig. 3. The limb profile of the hurricane Isabel area as measured by 
SCIAMACHY (the vertical resolution ~3 km). 


optical thickness is shown in Fig. 4(c). The most fre- 
quent values of t are around 5. The dependence as 
shown in Fig. 4(c) was also found analysing a great 
number of ground-based measurements by Trishchenko 
et al. (2001). This indirectly confirms the results ob- 
tained. SACURA does not retrieve values of t > 100 
because of insensitivity of the cloud reflection function 
to cloud optical thickness for extremely thick cloud sys- 
tems. So retrieved values of t larger than 100 are artifi- 
cially set to t= 100. It follows from Fig. 4(c) that it is 
about 50 such cases in the orbit analysed. Most of them 
are related to the hurricane Isabel. The liquid water path 
distribution obtained for the same orbit is given in Fig. 
4(d). Large values of LWP are again associated with the 
hurricane. However, we see that most frequent values of 
LWP are located in the region 50-70 g/m”. The average 
value of LWP is around 100 g/m?. This closely coincides 
with average annual values of LWP derived by Han 
et al. (1994) (87 g/m* over ocean) from AVHRR data 
on board NOAA-9 and NOAA-10. 

The measured top-of-atmosphere reflection function 
statistical distribution at the wavelength 440 nm is given 
in Fig. 5. It is obtained by the normalization of the nadir 
backscattered intensity by the solar light flux reflected 
by the white diffuse plate present in SCIAMACHY. 
The values of R in the range 0.3-1.0 are associated with 
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Fig. 4. The phase index (a), the cloud top height (b), the optical thickness (c), the liquid water path (d) statistical distributions for the orbit 08095 
(September 17th, 2003) as derived from SACURA using SCIAMACHY data. 


highly reflecting surfaces. So the study of the spatial dis- It follows that values of R below 0.1 never occur. This 
tribution of the reflection function is a convenient tool is because light reflectance from ocean/land, Rayleigh 
to access the reflective properties of observed surfaces scattering, and background aerosol are always present 


and atmospheric layers from space. in the satellite-detected signal at wavelengths used. 
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Fig. 5. The statistical distribution of the reflection function for the 
orbit 08095 (September 17th, 2003) as measured by SCIAMACHY 
(the spatial resolution ~30*60 km”). 


Three modes are clearly seen in Fig. 5. They can be asso- 
ciated with cloudless conditions (the mode for smallest 
values of R) and two modes for thin and thick (most 
probably, multilayered) clouds. It was found that the 
average value of R for the orbit studied is 0.45 with 
the coefficient of variance 50%. 


6. Conclusions 


The system of analytical equations has been devel- 
oped, which was used to retrieve the effective radius of 
droplets and the liquid water path of optically thick 
clouds (t > 5). Also other cloud characteristics are 
found. They are the optical thickness, the cloud albedo, 
the droplet column concentration, the phase index, the 
cloud top height and the cloud geometrical thickness. 

We give maps of the optical thickness, the effective ra- 
dius, the liquid water path, the phase index, R(443 nm), 
and cloud mask at our website for a particular orbit. 
These files are accompanied by the numerical data. 

The cloud top height is determined using spectral 
measurements in the 755-770 nm oxygen A-band and 


the phase index is determined using differences in the 
refractive index of liquid water as compared to ice in 
the spectral range 1550-1670 nm. Such spectral mea- 
surements are routinely performed by the SCIAM- 
ACHY instrument (Bovensmann et al., 1999). 

One important feature of the algorithm proposed is 
its simplicity and a high speed of retrieval procedures. 
It can be used as an integrated and autonomous part 
of operational cloud satellite retrieval schemes. It is 
characterized by a high accuracy and have important 
advantages as compared to the standard fitting methods 
used in up-to-date cloud retrieval schemes (especially if 
the speed of calculations is a concern as it is often the 
case in the analysis of large databases). The accuracy 
of the retrievals was checked by Kokhanovsky et al. 
(2004), Rozanov and Kokhanovsky (2004), and Kokha- 
novsky et al. (2005a,b). Nauss and Kokhanovsky (2005) 
found that SACURA produces results comparable to 
other algorithms (e.g., JAXA and NASA algorithms 
as applied to MODerate resolutlon Spectrometer 
(MODIS) data). 

Although both aerosol and molecular scattering and 
absorption are neglected, their influence is minimized 
by a careful selection of spectral channels. Also, the 
influence of molecular and aerosol scattering is of 
importance mostly for t<5 (Wang and King, 1997). 
This range of optical thicknesses is out of scope of the 
applicability of the algorithm described. The sensitivity 
of the top-of-atmosphere radiance to the ground albedo 
is low for thick clouds. So we assume that the albedo of 
ground is equal to zero. 

The applicability of the algorithm proposed for the 
optical thickness retrieval is limited both from the side 
of small (t<5) and large (t > 100) values of the 
optical thickness (Kokhanovsky et al., 2003). Note 
that the dependence of R on t disappears as t — oo. 
However, the appearance of such clouds in the Earth’s 
atmosphere is rare as indicated by Trishchenko et al. 
(2001). For the effective radius retrieval, however, 
there is no upper boundary for large t. This is due 
to the fact that the reflection of light depends on 
the size of droplets even for an infinitely thick 
absorbing cloud. The footprint of the SCIAMACHY 
instrument is 30*60 km? for all studies performed 
here. This underlines the fact that thin clouds retri- 
evals with this instrument is difficult to handle due 
to possible broken cloud conditions associated with 
small values of t. 

In conclusion, note that the colour index c can be 
used for mapping convective cloud events using 
SCIAMACHY limb measurements (see Fig. 3). This is 
of importance for a number of problems and especially 
for troposphere-stratosphere exchange studies. 

SACURA code is available online as a part of the 
radiative transfer solver SCIATRAN (Rozanov et al., 
2005) (see www.iup.physik.uni-bremen.de/sciatran). 
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